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We study entanglement witnesses that can be constructed with regard to the 
geometrical structure of the Hilbert-Schmidt space, i.e. we present how to use these 
witnesses in the context of quantifying entanglement and the detection of bound 

<^ • entangled states. We give examples for a particular three-parameter family of states 

O . that are part of the magic simplex of two-qutrit states. 

The determination whether a general quantum state is entangled or not is of utmost 
importance in quantum information. States of composite systems can be classified in Hilbert 
space as either entangled or separable. The geometric structure of entanglement is of high 
interest, particularly in higher dimensions. For 2x2 bipartite qubits, the geometry is 
highly symmetric and very well known. In higher dimensions, like 3x3 two-qutrit states, 
the structure of the corresponding Hilbert space is more complicated and less known. New 
phenomena like bound entanglement occur. In principle, all entangled states can be detected 
by the entanglement witness procedure. We construct in this Brief Report entanglement 
witnesses with regard to the geometric structure of the Hilbert-Schmidt space and explore 
a specific class of states, the three-parameter family of states. We discover by our method 
^ \ new regions of bound entanglement. 

t^j- ■ Let us recall some basic definitions we need in our discussion. We consider a Hilbert- 

Schmidt space A = Ai <E> A% ® ... ® A n of operators acting on the Hilbert space of n 
composite quantum systems - it is of dimension d\ x <i 2 x • • • x d n and D := d\d 2 . . . d n . 
In Ref. [H a new geometric entanglement measure (which we call Hilbert-Schmidt measure) 
based on the Hilbert-Schmidt distance between density matrices is presented - and discussed 
in Ref. 0] - that is an instance of a distance measure (see Refs. @, 0|). The entanglement 
of a state p can be quantified (or "measured") via the minimal Hilbert-Schmidt distance of 
the state to the convex and compact set of separable (disentangled) states S: 



o 



D(p) :— min ||cr — p|| . (1) 



It is defined on the Hilbert-Schmidt metric with a scalar product between operators that 
are elements of the Hilbert-Schmidt space A: (A, B) := TrA^B, A, B G A and the norm 
\\A\\ := a/ (A, A) . Of course these definitions apply to density matrices (states of the 
quantum system), since they are operators of A with the properties pt = p, Tip = 1 and 
p > (positive semidefinite operators). 

Because the norm is continuous and the set of separable states S is compact the minimum 
in Eq. ([TJ is attained for some separable state <t , min^s ||a — p\\ = \\a — p\\ , which we 
call the nearest separable state to p. Clearly if p G S then a = p and D(p) = 0. If p is 
entangled, Pen^ then D(p) > and cr lies on the boundary of the set S. It is shown in 
Refs. p, [y, [71, Isj] that an operator 



Aipt — CTq — p cnt — ((To, CT — p cn t) Id 



(2) 
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defines a hyperplane including the state a that is tangent to S, and it has the properties 

(p cnt , A opt ) = Tr p cnt A pt < , (3) 
(a,A opt ) = Tr crA opt > Wa E S . (4) 

It is clearly Hermitian and thus is an entanglement witness for the state p e nt @,[Iq|. More- 
over, since the hyperplane is tangent to S (Tro"oA op t = 0) the operator A opt (j2J) is called 
optimal entanglement witness. We call any entanglement witness that is constructed in the 
way of Eq. (j2J) - with any two states, not necessarily the nearest separable state - a geometric 
entanglement witness. 

The crucial point lies in finding the nearest separable state: Once found, we can both 
quantify the entanglement of a state and construct an (even optimal) entanglement witness. 
But finding the nearest separable state is a hard task. A "guess-method" , that is a method 
to check if a good guess for the nearest separable state is indeed right, is presented in Ref. jl]. 
Another way is to look for the nearest state that is positive under partial transposition (PPT) 
first. A method that finds the nearest PPT state in many cases is presented in Ref. 

The PPT-criterion j9|, [l2| is a necessary criterion for separability (sufficient for 2 x 2 or 
2x3 dimensional Hilbert spaces): A separable state has to stay positive semidefinite under 
partial transposition - i.e. transposition in only one subsystem. Thus if a density matrix 
becomes indefinite under partial transposition, i.e. one or more eigenvalues are negative, it 
has to be entangled and we call it a NPT entangled state. But there exist entangled states 
that remain positive semidefinite - PPT entangled states - these are called bound entangled 



states, since they cannot be distilled to a maximally entangled state [131 . 114 . 

The set of all PPT states P is convex and compact and contains the set of separable 
states. Thus the nearest separable state cr can be replaced by the nearest PPT state To for 
which the minimal distance to the set of PPT states is attained, min re p \\r — p\\ = ||r — p|| . 
If p is a NPT entangled state Pnpt and r the nearest PPT state, then the operator 

App T := r — Pnpt — (to,t — Pnpt)1d (5) 

defines a tangent hyperplane to the set P for the same geometric reasons as operator (J2J) 
and has to be an entanglement witness since P D S. In principle the entanglement of 
Pnpt can be measured in experiments that should verify TrAp PT p NPT < 0. If the state r 
is separable, it has to be the nearest separable state a since the operator (jSJ) defines a 
tangent hyperplane to the set of separable states. Therefore in this case App^ is an optimal 
entanglement witness, App^ = A opt , and the Hilbert-Schmidt measure of entanglement can 
be readily obtained. If tq is not separable, that is PPT and entangled, it has to be a bound 
entangled state. 

Unfortunately it is not trivial to check if the state To is separable or not. It is hard to 
find evidence of separability, but it might be easier to reveal bound entanglement, not only 
for the state To but for a whole family of states. A method to detect bound entangled states 
we are going to present. 

Consider any PPT state pppx and the family of states p\ that lie on the line between 
Pppt and the maximally mixed (and of course separable) state ^Idj 

p A := A p PPT H — — 1 D . (6) 

We can construct an operator C\ in the following way: 

C\ = Pa - Pppt - (px, Px - Pppt) Id • (7) 
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FIG. 1: Sketch of the presented method to detect bound entanglement with the geometric entangle- 
ment witness C\. The dashed line indicates the detected bound entangled states p\ (A m i n < A < 1), 
states p\ with < A < A m ; n can be separable or bound entangled (straight line). 



If we can show that for some A m ; n < 1 we have Tr aC\ . > for all a G S, C\ . is an 

mill ^min — " ^mm 

entanglement witness and therefore pppx and all states p\ with A m i n < A < 1 are bound 
entangled (see Fig. 

As an example we introduce the following family of three-parameter two-qutrit states 
(dimension 3x3 and 1 := 1 9 ): 

1 — Ol — 3 — T 3 T 

PcPa ■= -1 + «P o + | (Pio + P20) + j (P01 + Pu + P21) , (8) 

where the parameters are constrained by the positivity requirement P gjL- y > 0. The states 



P nm are projectors onto maximally entangled two-qutrit vector states |15|] - the Bell states 



Pnm ■= (U nm <g> 1)|0+) (<f>t\ (Ul m ® 1) , (9) 

where denotes the maximally entangled state = ^= J2j \ j) ® |j) an d t^nm represent 
the VFeyi operators U nm = X]fc=o e ^ kn \ k) ((k + m ) modd| with n, m — 0, 1, . . . , d — 1 (here 

^1=^2 = ^)- 

The states (jHJ) lie in the magic simplex of two-qutrit states which is the set of all two-qutrit 
states that can be written as a convex combination of the projectors P nm (Q. Viewing the 
indices nm as points in a discrete phase space the ma gic s im p lex reveals a high symmetry 



that makes its geometry much more evident (see Refs. m Il7|,ll8l). 

For 7 = the geometry of the states ([8]) becomes rather simple. It is shown in Ref. [If 
that all states of this two-parameter family are either NPT entangled or separable, that 
is, all PPT states coincide with the separable states. Therefore we can easily calculate the 
Hilbert-Schmidt measure with help of our entanglement witness (J2]). 

In the two-dimensional picture (see Fig. [2]) we can easily acquire the nearest points on the 
border of the separable states - here equal to the PPT states - to the points in the regions 
I and II that correspond to the entangled states: For region I the nearest point to (a, 3) is 
(4 + 0) which characterizes a state dp of the family (jHJ). For region II the nearest point is 
(^(—2 + 20a + 53), |(2 + 4a + 3)) which corresponds to a state a a ^. According to Ref. [s| 
we can check if a proposed separable state a is the nearest separable one to an entangled 
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±(P,o + P 2 o) 



FIG. 2: Illustration of states p Qj/ g © (7 = 0) and their partial transposition. The regions I and II 
label the regions where the states are NPT entangled, they are PPT and separable in the overlap 
with the region of PPT points. The PPT points become semipositive under partial transposition. 



state p (considering the geometry of all states) by checking if the operator 

C = a-p-(a,a-p)l B (10) 
is an entanglement witness. To do so let us first state the following Lemma: 

Lemma 1. For any Hermitian operator C of a bipartite Hilbert- Schmidt space of dimension 
d x d that is of the form 



d-l 



C = a { (d - 1) l d 2 + ^2 c «™ u nm ® U„ nm , a e R + , c nm e C 



n,m=0 



the expectation value for all separable states is positive, 

(p, C) > V p G S , if \c nm \ < 1 Vn,m. 
Proof. Any bipartite separable state can be decomposed into Weyl operators as 



(11) 



(12) 



^ d-l d-l 

J2,Pk ^(ld ® Id + Yl ^ d ~ 1 U nm Unm ® Id + ^ Vd - 1 m? fc l d ® I7 tt 



n,m=0 



Lk=0 



di-1 



+ ( d - 1 ) 7l nm^ffc^«m®^ 
n,m,i,fc=0 



/A 



|n fc | < 1 , \m k \ < 1 



>0, $> fc = l 



(13) 
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where we define |ra fe | 2 := J2nm n nm n nm- We call a decomposition (TT5j) Block vector form of 
the density matrix. Note that for d > 2 not any vector with complex components represents 
a Bloch vector (a quantum state), details can be found in Ref. [191 ]. 

Performing the trace we obtain (keeping notation p< formula (JHj) becomes more evident) 



(p,C) = TrptC = 2p fc f(d-l)afl+^c 
and using the restriction |c nm | < 1 Vn, m we have 



(14) 



L -nm"'nm"''-nm 



< J2\ n nmW m -nm\ < 1, (15) 



and since the convex sum of positive terms stays positive we get (p,C) > Vp G S . □ 

We set up the operators Cj, C/j for the two regions, respectively, according to Eq. (JTUJ) 
and write them in terms of Weyl operators normalized by ||<7 — p\\ for convenience, 

C I = -^=(2t-U 1 -U 2 ), C IJ = -±={21 + U l -U 2 ), (16) 



where 



Ui := U 01 <g> C/oi + C^02 ® ^02 + C/11 ® E/Lu + U 12 g> C/_ia + E/21 ® J7_ 2 i + ^22 ® E 7 - 
C/ 2 := ^2 + with f/ 2 7 := C/ 10 ® E7_ 10 , f/ 2 7/ := C/20 ® ^-20 • 
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(17) 



Both operators Cj and Cu satisfy Eq. (1121) of Lemma [TJ Due to the construction of the op- 
erators, Eq. ([3]) is satisfied for any entangled states p a ,p in the regions I and II. Consequently 
Ci and Cu are entanglement witnesses and therefore dp and a a ,p are the nearest separa- 
ble states Co;/? and cro- a ,p- The corresponding Hilbert-Schmidt measures of the entangled 
two-parameter states p a% p are 



2^2 ( 1 1 



^/(PS) = IK/» - PSII = ^ ( « - 7 - ) • (18) 



2^/2 ( 1 5 



Note that the measures can also be viewed as a maximal violation of the entanglement 
witness inequality PJ, as it is shown in detail in Refs. @, H|. 

Another way to arrive at the nearest separable states for the two-parameter states is to 



calculate the nearest PPT states with the method of Ref. 11 first and then check if the 



gained states are separable. If we do so we obtain for the nearest PPT states the states dp 



and a a ,p we have found with our "guess" method, we know from Ref. [16J] that these states 
are separable and therefore they have to be the nearest separable states. 

Let us return to the family of three-parameter states p a ,p,-y ©• For 7 7^ it is not trivial 
to find the nearest separable states since the PPT states do not necessarily coincide with the 
separable states. But we can use our geometric entanglement witness ([7]) to detect bound 
entanglement. 
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FIG. 3: The states Pafi,^ are depicted with slices of the three-dimensional parameter space for 
fixed values of 7. The curved borders of the PPT regions for 7^0 are hyperbolae that result from 
intersecting a cone of PPT points with the planes of constant 7. The slices for positive and negative 
7 overlap. All parameter axes are chosen non-orthogonal such that they become orthogonal to the 
boundary of the positivity region in order to reproduce the symmetry of the magic simplex. 



In Ref. [20j the following one-parameter family of two-qutrit states was introduced: 

Pb = ~|<^)(<^| + ~a + + ^ ( 7_, 0<6<5, (20) 

where 

o+ := 1/3 (|01) (01| + |12) (12| + |20) (20|) , 

ct_ := 1/3 (|10) (10| + |21) (21 1 + |02) (02|) . (21) 

Let us call this family of states Horodecki states. Interestingly, the states (|2"U|) are part of the 
three-parameter family (|SJ), namely p^ = p a ,0,-y > with a — (6 — b)/21 , (3 = —26/21 , 7 = 
(5 — 26) /7 and thus lie in the magic simplex. Testing the partial transposition we find that 
the Horodecki states ([20]) are NPT for < b < 1 and 4 < b < 5 and PPT for 1 < b < 4. 
In Ref. [2(| it is shown that the states are separable for 2 < b < 3 and bound entangled for 
3 < b < 4. 

We now want to pursue the following idea: Starting from a PPT and entangled - bound 
entangled - Horodecki state pf we construct operators in the way of Eq. ([7]) and try to find 
bound entanglement on the line between pf E and the maximally mixed state - see Fig. [TJ 

For a geometric picture of the three-parameter family of states p a ^ n we can fix values 
of 7 and draw two-dimensional slices (see Fig. [3]). Translated to values of 7 the states 
Pb = Pa,/3,-y are bound entangled for —3/7 < 7 < — 1/7 and in each slice of this region 
there is exactly one point that corresponds to the bound entangled Horodecki state, namely 
(cv = (1 + 7)/6 , (3 — (—5 + 77)/21). The one-parameter line of the states pb therefore cuts 
through the slices of fixed 7 and lies on the boundary of the states ([8]). 
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For a fixed 7 the operator C 7i a © has the following form: 

C 7 ,a = px-pf E -(px, P x-pf E )l = a(2t + c l U 1 + c 2 U^ + 4Ui I ), 

l + 37 2 w , 1N 8 2(1-7^7*) 

with 0= -A(A-l), c x = rr-, c 2 = ^44, 22 

36 v ; ' 7A(l + 3 7 2 )' 7A(l + 3 7 2 ) ' v ' 

and p\ represents the line of states 

Px = \fi B + ^l. (23) 

It is convenient to define functions /i(7, A) := |ci| and ^2(7, A) := \c 2 \ ■ We are interested to 
find the minimal A for a particular 7 which is attained at max{/i(7, A), ^(7, A)} = 1, i.e. the 
minimal A such that C 7t \ f !22p is necessarily an entanglement witness (see Lemma [TJ). There 
exist values of A with A < 1 such that max{/i(7, A), f 2 (j, A)} < 1 for l/y/21 < \^y\ < 3/7. 
That means we are able to detect bound entangled states on lines ( 123!) between the bound 
entangled Horodecki states and the maximally mixed state |l until a value of I7I = l/y/21. 
These lines form a planar section of bound entangled states restricted by A m i n < A < 1. 
As a "side product" we are able to detect bound entanglement for the Horodecki states for 
1 < b < |(15 - v / 2T)(^ 1-74) and ±(15 + y/21){~ 3.26) < b < 4. 

A detailed examination of the coefficient functions /i(7, A) and ^(7, A) exhibits that 
for l/y/21 < I7I < a/5/ 7 we have fx > f 2 and for a/5/7 < I7I < 3/7 we have fx < f 2 and 
equality fx = f 2 at v5/7. The total minimum is A^* n = 7/8 = 0.875, i.e. the minimal A 
that can be reached such that all states p\ (j23"j) are detected to be bound entangled for 
A£* n < A < 1. It is attained for | 7 1 = a/5/7 or b= |(5-a/5) - 1.38 and 6 = |(5+a/5) ~ 3.62. 

Summarizing, the magic simplex contains three classes of states: NPT entangled states, 
separable states but also PPT entangled or bound entangled states. We have found such 
bound entangled states analytically in a quite large region of the magic simplex. Our method 
works very well for the examples presented and is geometrically very intuitive. However, we 
cannot detect in these examples the border between bound entangled and separable states. 
The reason is that condition \c nm \ < 1 Vn,m in Eq. (|T2|) is only sufficient for C being an 
entanglement witness and not necessary. There is generally still no operational method for 
detecting the nearest separable states. 
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